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Abstract 

In this work, we present the gravitational field generated by a cos- 
mic string carrying a timelike current in the scalar-tensor gravities. 
The mechanism of formation and evolution of wakes is fully investi- 
gated in this framework. We show explicitly that the inclusion of elec- 
tromagnetic properties for the string induces logarithmic divergences 
in the accretion problem. 
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1 Introduction: 



Topological defects are predicted in many gauge models as solitonic solutions 
resulting from spontaneous breaking of gauge or global symmetries. Among 
all these solutions, cosmic strings have attracted attention because they may 
be the source of large-scale structure in the Universe [1]. A relevant mecha- 
nism to understand the structure formation by cosmic strings involves long 
strings moving with relativistic speed in the normal plane, giving rise to ve- 
locity perturbations in their wake [2]. The mechanism of forming wakes has 
been considered by many authors in the General Relativity theory [3, 4]. 
More recently, Masalskiene and Guimaraes [5] have considered the formation 
of wakes by long strings in the contexte of scalar-tensor theories. They have 
shown that the presence of a gravitational scalar field - which from now on we 
will call generically as "dilaton" - induces a very similar structure as in the 
case of a wiggly cosmic string [4]. Further, Bezerra and Ferreira [6] showed 
that, if a torsion is presented, this effect is amplified. 

Our purpose in this paper is twofold. We first generalize the works of 
Peter and Puy [7] and Ferreira et al. [8] for the case of a cosmic string carrying 
a timelike current in the framework of scalar-tensor gravities. From the point 
of view of purely gravitational physics, it was shown that the inclusion of a 
current hardly affects the metric outside the string [9]. However, it is well- 
known that inclusion of such an internal structure can change the predictions 
of cosmic string models in the microwave background anisotropies [10, 11]. 
This is our second goal. Namely, we analyse the formation and evolution of 
wakes in this spacetime and we show explicitly how the current affects this 
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mechanism. For this purpose, we consider a model in which non-baryonic cold 
dark matter propagates around the current- carrying string. The Zel'dovich 
approximation is carried out in order to treat this motion. We anticipate that 
our main result is to show that the inclusion of a current brings logarithmic 
divergences and can actually break down the accretion mechanism by wakes. 

This work is outlined as follows. In section 2, after setting the rele- 
vant microscopic model which describes a superconducting string carrying a 
current of timelike type, we present its gravitational field and we show the 
explicit dependence of the deficit angle, the deflection of light and geodesies 
on the timelike current. In section 3, we consider the mechanism of forma- 
tion and evolution of wakes in this framework by means of the Zel'dovich 
approximation. Finally, in section 4, we end up with some conclusions and 
remarks. 

2 Timelike Current-Carrying Strings: 

In this section we will study the gravitational field generated by a string 
carrying a current of timelike-type. We start with the action in the Jordan- 
Fierz frame 

s = itb/ d4x ^m-^r u 9^dM+ s m [^~g,u] (i) 

g^v is the physical metric which contains both scalar and tensor degrees of 
freedom, R is the curvature scalar associated to it and S m is the action 
for general matter fields which, at this point, is left arbitrary. The metric 
signature is assumed to be (— , +, +, +). 

3 



In what follows, we will concentrate our attention to superconducting 
vortex configurations which arise from the spontaneous breaking of the sym- 
metry U(l) x U em (l). Therefore, the action for the matter fields will be 
composed by two pairs of coupled complex scalar and gauge fields (ip, B^) 
and (a, A^). Also, for technical purposes, it is preferable to work in the so- 
called Einstein (or conformal) frame, in which the scalar and tensor degrees 
of freedom do not mix. 

S = lQ^J d ' Xv ^ [R ~ 2d ^ d " <P] 

+ J d'x^-gl-^miD^TD^ + {D^fD^a) 

- ^F, u F^ + H, u Hn-^ 2 ^)V{MM\)], (2) 

where = dfiA u — d u A^ , H^ v = djj,B v — d u B^ and the potential is suitably 
chosen in order that the pair (</?, B^) breaks one symmetry U{\) in vacuum 
(giving rise to the vortex configuration) and the second pair (a, A^) breaks the 
symmetry U em {l) in the core of the vortex (giving rise to the superconducting 
properties) 

v(\<p\, k|) = ^(M 2 - v 2 ) 2 + /(M 2 - v 2 )W\ 2 + ^H 4 + ^H 2 (3) 

We restrict, then, ourselves to the configurations corresponding to an 
isolated, static current-carrying vortex lying in the z — axis. In a cylindrical 
coordinate system (t, r, z, 6) such that r > and < 6 < 2%, we make the 
following ansatz: 

V = v(r)e w Bfl = -[Q(r)-l] (4) 
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The pair (a, A^), which is responsible for the superconducting properties of 
the vortex, is set in the form 

a = a(r)j*® A t = -[P t (r) - d tX ] (5) 

e 

where Pt corresponds to the electric field which leads to a timelike current 
in the vortex. We also require that the functions (p, Q(r), <r(r)andPi must be 
regular everywhere and must satisfy the usual boundary conditions of vortex 
[12] and superconducting configurations [13, 14]. 

The action (2) is obtained from (1) by a conformal transformation 

~g ilv = n 2 ((j))g^ 

and by the redefinition of the quantity 

which makes evident the feature that any gravitational phenomena will be 
affected by the variation of the gravitation constant in the scalar-tensor 
gravity, and by introducing a new parameter 

Variation of the action (2) with respect to the metric g^ v and to the 
dilaton field <fi gives the modified Einstein's equations and a wave equation 
for the dilaton, respectively. Namely, 

U g( j) = -4nG.a(<f>)T (6) 
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Where T^v is the energy-momentum tensor which is obtained by 



T = 



-2 6S, 



mat 



(7) 



V~9 &9ixv 

We note, in passing, that, in the conformal frame, this tensor is not conserved 
providing us with an additional equation 

V M T^ = a{4>)TVA- 

In what follows, we will write the general static metric with cylindrical 
symmetry corresponding to the electric case in the form 



(8) 



where ip, 7, (3 are functions of r only. 

The non-vanishing components of the energy-momentum tensor using (7) 

are 

T* = -ifi 2 (0) {-g*o>P t 2 + <T [p' 2 + a' 2 ] + g ed v 2 Q 2 } 



P 

At 1 t 



a 



8tT 



+ 



fi Jn2 rr 



-£l\^)V(a, V ) 



t; = ~^){g tt ^P t 2 -g 



if + <T 



,2 }+g e \ 2 Q 2 } 



1 f P' 2 O' 2 



8tt' 



1 

_< 

2 

1 

8tt' 



l g = {<p)[g a P t +g |y? + a j - g ip Q j 

e 2 



(9) 



^ = -^ 2 (0){^VP, 2 + ^[/ + a' 2 ]+/VQ 2 } 



1 [ P' z 



-fi 4 (0)V>,^). 
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Therefore, for the electric case, eqs. (6) are written as 

(PiPj = 4vrG*e 2(7 ^ ) / 3[-T + + T 2 + T 3 3 ] 

(Ptj/y = -AnG^e^-^T (10) 

In order to solve the above equations we will divide the space in two 
regions: the exterior region, r > r , in which only the electric component 
of the Maxwell tensor contributes to the energy-momentum tensor and the 
internal region, < r < r , where all matter fields survive. r is the string 
thickness. 

2.1 The Exterior Metric: 

Due to the specific properties of the Maxwell tensor 

T; = and T^ = ^(T^T afi )SS (11) 

the Einstein's equations (10) may be transformed into some algebraic rela- 
tions called Rainich algebra [15, 16] which, for the electric case, have the 
form 1 

R = R z z + R r r + R e e + R t t = 2e 2{ ^ ) <P' 2 

[R\f = (R r r - 2g rr <f>' 2 ) 2 = {R e e ) 2 = (R z z ) 2 

R\ = — Rq R^ = R r r — 2g rr <\j R® = R z z (12) 
^dn the scalar-tensor theories, these relations are modified by a term which depend on 
the dilaton [8]. 
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Then we have 



the solution is 



+ V-^ ,2 -7" = ^ (13) 

7" + -i = (14) 
r 



(7* 
— 

* = " ln (9- ln (trr) (15) 



Therefore, the exterior metric for a timelike current-carrying string is: 

( r \ 2m 2 — 2n / r \ —2n 

- ) W\r){dr 2 + dz 2 ) + ( - ) W 2 (r)B 2 rW 

To J V r / 



r / Vr , 

2n -| 

rft 2 (16) 



r J W 2 (r) 



where 



W{r) ee V "/ +fc . (17) 



The constants m, n, k, D will be determined after the inclusion of the matter 
fields. 



2.2 The Internal Metric and the Matching: 

In order to solve the Einstein's equations (6) in the internal region (where 
all matter fields will contribute to the energy-momentum tensor) we will use 
the weak-field approximation. Therefore, we can expand the metric and the 
dilaton fields to first order in G such that 

(f) = o + 0i 
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In this way, eqs. (6) reduce to 

V 2 V = -16ttGo(tW - \n, v T^). (18) 

in a harmonic coordinate system such that (h% — \5^h) >v = 0, and the lin- 
earised equation for the dilaton field 

V 2 0(i) = -4^G o a(0 o )T (o) . (19) 

In the above eqs. (18) and (19), TjfJ is the energy-momentum tensor eval- 
uated to zeroth-order in Go, and its trace. In cartesian (harmonic) 
coodinates the nonvanishing componentes of Tf£) are 

rpt ^ / Jn2 | , /2 , /2 1 ,„2 ri 2\ 

T (o)t = { a p t + V + a + ^ Q ) 
8tt { e 2 ^ r 2 g 2 ~*~} ( ' ^) 



I5,„ = ( S in 2 (fl)-i)L' 2 + </ 2 



t 2 o 2 p; 2 



r 2 e47r 2 

= -l{-* 2 P t 2 + <p' 2 + ° ,2 + ^ 2 Q 2 } 

"' +^\-V(<r 1 ip) (20) 



8n \ e 2 r 2 q 2 



In this paper, in order to solve eqs. (18) and (19) with source given by (20), 
we will use a method which was first applied by Linet in ref. [17] and later 



by Peter and Puy [7] and Ferreira et al. [8]. The trick consists in writing 
down the string and its components fields by means of S functions 2 . In doing 
this, we are allowed to write the string's energy-momentum tensor as 

T» v = diag (U, -Z, -Z, —T) (21) 

where U, Z and T are macroscopic quantities which are defined as the energy 
per unit length, the transverse components per unit length and the tension 
per unit length, respectively: 

U = Mj = -2tt I" Tl r dr 



t 

ro 



rro 

X = M r r = -2tt / T r r r dr (22) 
Jo 

Y = M e d = -2ir / T e d r dr 
Jo 



rro 

T = Ml — -2tt / T" r dr 



rr 

-27T / 

Jo 



From the transverse quantities X and Y we define one single cartesian com- 
ponent Z: 

Z=- J rdrdQTl = - J rdrddT y y (23) 

On the other hand, the energy-momentum in the exterior region (and 
corresponding definition of the electric current) can be re-written in cartesian 
coordinates as 

I 2 



rptt rpzz 

em - em - ^ 



T e k l = -^- 4 (2x k xi-r 2 5 kj ) (24) 



2 The calculations are straightforward but lengthy. We will skip here the details and 
provide directly the results. For the details of these calculations, we refer the reader to 
refs. [7, 8] . 
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Therefore, we have 



rptt 

I (0) 



rpZZ 

I (0) 



rpkj _ 

i (o) - 



US(x)S(y) + — V 2 
-T5(x)5(y) + 



S kj S{x)5{y) 



In 



r 



In I - 



1 2 



^0 

d k d 3 In (^) 



(25) 



2tt 



Solving the linearised Einstein's equations, we obtain 



V 2 h tt = -4G 
K = -4G 



r 



I 2 V 2 In - + 27r(C/ — T — I )5(x)5(y) 



12 HQ) ^ u - t -* 2 H 



(26) 



V 2 ^ = -4G 
/i** = -4G 



I 2 V 2 (in ^ ) 2 + 2tt(^ - T + I 2 )5(x)5(y) 
I 2 (^))\ i U-T + I 2 ) ^ 



(27) 



V 2 /i fcj 



-4G 



+2/ 2 <9 fc <9 i In (-^) + 27r^(C/ — T — I 2 )5(x)5(y) 



h kj = 2G I 2 r 2 d k d J In (J^j - AG S kj (U + T-I 2 ) In 



(28) 



^o/ Vr , 
In order to match the internal and the external solutions and obtain the 
parameters m, B, n e D we need to return to the cylindrical coordinates 
system 



9u 

9zz 



- 1 + 4G 



1-AG 



<^) 2 + {U - T - l2) <)}}> 
/ 2 (ln(^)) 2 + ^-T + / 2 )ln(- 



(29) 



(30) 
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g rr = 1 - 2G I 2 -4G (U + T- I 2 ) In 



1 + 2G r - AG (U + T - I 2 ) In 



/ r 



Vr 



(31) 
(32) 



Before presenting the final result we note that the metric (29) does not pre- 
serve the feature gnu = g zz . We need again to make a change of coordinates, 
to first order in Go 

l + ai -a 2 ln(Q-a 3 (inQ-)) 2 

where a x = G (4T - I 2 ), «2 = 4G T and o 3 = -2G I 2 . Without loss of 
generalities, let us make R = r and expand the parameters in powers of 
In (^)- In doing that, we finally get 



R = r 



ds 2 



1-AGo 



(U-T + I 2 ) lnQ-) + / 2 ln 2 ^)]}(dr 2 + dz 2 ) 



+ r l 



l-8G [T--j - 4G (C/ - T - I 2 ) ln(-J - 4G / 2 In 2 



r 



{l + 4G„[i» m>(Q + (U-T-l>) 



In the next section, we will use this metric to compute the deficit angle 
and the deflection of light rays by a current-carrying string. 



2.3 Deficit Angle and the Deflection of Light Rays by 
a Current- Carrying String: 

Now that we have calculated the metric (33) of a timelike current-carrying 
string, we can compute some quantities associated to it, such as the deficit 
angle and the deflection of light rays. The deficit angle is calculated with the 
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help of the expression 



5 e = 2tt(1 



1 d 



For the metric (33) it is, to first order in Gq: 



5 e = 4ttG (U + T- 2J 2 ) + O(G 2 ) 



(34) 



It is clear that, if I is zero (the neutral string case), the deficit angle (34) 
reduces to 8$ = 8ttG U which is of order ~ 10~ 6 for GUT strings. 
The deflection of the light ray is given by the expression 



From expression (36) we can observe some remarkable features of the de- 
flection of light due to a timelike current-carrying string in the scalar-tensor 
gravity. First of all, there is no dependence on the tension in the deflection of 
light rays. This feature is a reminescence of the neutral string. The second 
feature is the dependence in the impact parameter p due to the presence 
of the current. This is due to the fact that the inclusion of an electromag- 
netic field actually increases the total energy with the distance to the string 
core and, as a consequence, the deflection of light is more intense. Finally, 
we notice that the dilaton field amplifies this effect since its inclusion also 
contributes to increase the total energy of the system. 




(35) 



As in [7] we obtain, knowing that p is the impact parameter, 
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3 Formation and Evolution of the Wakes and 
the Zel'dovich Approximation: 



A relevant mechanism to understand the structure formation by cosmic 
strings involves long strings moving with relativistic speed in the normal 
plane, giving rise to velocity perturbations in their wake [2] . Matter through 
which a long string moves acquires a boost in the direction of the surface 
swept out by the string. Matter moves toward this surface by gravitational 
attraction and, as a consequence, a wake is formed behind the string. In what 
follows, we will study the implications of a timelike current on the formation 
and evolution of a wake behind the string which generates the metric (33). 
For this purpose, we will mimic this situation with a simple model consist- 
ing of cold dark matter composed by non-relativistic collisionless particles 
moving past a long string. In order to make a quantitative description of ac- 
cretion onto wakes, we will use the Zel'dovich approximation, which consists 
in considering the Newtonian accretion problem in an expanding universe 
using the method of linear perturbations. 

To start with, we first compute the velocity perturbation of massive par- 
ticles moving past the string. If we consider that the string is moving with 
normal velocity v s through matter, the velocity perturbation can be calcu- 
lated with the help of the gravitational force due to metric (33): 



u = %nGoUv s r ) 



(37) 



2a(0 o ) 2 (U + T + I 2 ) - (U - T - I 2 ) - 2 In [pj I 2 
with 7 = (1 — v 2 )^ 1 / 2 . The first term in (37) is equivalent to the relative 
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velocity of particles flowing past a string in general relativity. The other 
terms come as a consequence of the scalar-tensor coupling of the gravitational 
interaction and the superconducting properties of the string. 

Let us suppose now that the wake was formed at tt > t eq . The physical 
trajectory of a dark particle can be written as 

h(x, t) = a(t) [x + i/j(x, t)] (38) 

where x is the unperturbed comoving position of the particle and ip(x,t) is 
the comoving displacement developed as a consequence of the gravitational 
attraction induced by the wake on the particle. Suppose, for simplification, 
that the wake is perpendicular to the x-axis (assuming that dz = in the 
metric (33) and r = \Jx 2 + y 2 ) in such a way that the only non-vanishing 
component of ip is ip x . Therefore, the equation of motion for a dark particle 
in the Newtonian limit is 

h = -V fc $ (39) 
where the Newtonian potential $ satisfies the Poisson equation 

= AtvG oP (40) 

where pit) is the dark matter density in a cold dark matter universe. For 
a flat universe in the matter-dominated era, a(t) ~ t 2 ^ 3 . Therefore, the 
linearised equation for if) x is 

^ + 7^-J^ = (41) 
with appropriated initial conditions: ip(U) = and = —Ui- Eq. (41) is 
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the Euler equation whose solution is easily found 



i(>(x,t) = - 
o 



1 , t ^ 2/S 



Calculating the comoving coordinate x(t) using the fact that h — in the 
"turn around" 3 , we get 



x(t) 



Uitf 



Uiti — 



t ,2/3 



(42) 



t Vti 

With the help of (42) we can compute both the thickness d(t) and the surface 
density a(t) of the wake [1]. We have, then, respectively (to first order in 
G ) 



d(t) 



— - MG Uv sl + — °- 
5 \UJ { v s -f 

2 i /n 1/3 r 



2a((p ) 2 (U + T + I 2 ) - (U - T - I 2 ) - 2 In (-^ / 



(^) [S^(^ 7 ) 2 + 2«(0 O ) 2 (U + T + / 2 ) - (u - T - I 2 ) - 2 In (^-) / 



(Tit) W ( - ) 



where we have used the fact that p{t) 



(43) 



i 



6nG t 2 



for a flat universe in the 



matter-dominated era and that the wake was formed at U ~ t eq . Clearly, 
from eq. (43) we see that the presence of the current makes the accretion 
mechanism by wakes to diverge. 



4 Conclusion: 

Inclusion of a current in the internal structure of a cosmic string could dras- 
tically change the predictions of such models in the microwave background 
3 The moment when the dark particle stops expanding with the Hubble flow and starts 
to collapse onto the wake. 
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anisotropies. In particular, a current of a timelike-type could bring some 
divergences leading to some unbounded gravitational effects [9] . In this work 
we studied the effects of a timelike-current string in the mechanism of wakes 
formation. For this purpose, we first studied the gravitational properties of 
the spacetime generated by this string in the framework of scalar-tensor grav- 
ities. We analysed the dependence of the deficit angle and the deflection of 
light on the current in this spacetime. Then, we carried out an investigation 
of the mechanism of formation and evolution of wakes in this framework, 
showing the explicit contribution of the current to this effect. 

Wakes produced by moving strings can provide an explanation for fila- 
mentary and sheetlike structures observed in the universe. A wake produced 
by the string in one Hubble time has the shape of a strip of width ~ v s ti. 
With the help of the surface density (43) we can compute the wake's linear 
mass density, say fi, 

<r(t) « ~ (£) V3 [8U(v sl f + 2a(0 o ) 2 (U + T + J 2 ) - (u - T - I 2 ) - 2 In 

(44) 

If the string moves slower or if we extrapolate our results to earlier epochs, 
we see that the logarithmic term would bring divergences and the mechanism 
of forming wakes would break down. 
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